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Abstract.  A  scheme  for  introducing  mountain  wave-induced 
temperature  pertubations  in  a  microphysical  PSC  model  has 
been  developed.  A  data  set  of  temperature  fluctuations  at¬ 
tributable  to  mountain  waves  as  computed  by  the  Moun¬ 
tain  Wave  Lorecast  Model  (MWLM-2)  has  been  used  for 
the  study.  The  PSC  model  has  variable  microphysics,  en¬ 
abling  different  nucleation  mechanisms  for  nitric  acid  trihy¬ 
drate,  NAT,  to  be  employed.  In  particular,  the  difference 
between  the  formation  of  NAT  and  ice  particles  in  a  sce¬ 
nario  where  NAT  formation  is  not  dependent  on  preexisting 
ice  particles,  allowing  NAT  to  form  at  temperatures  above 
the  ice  frost  point,  Tjce ,  and  a  scenario,  where  NAT  nucle¬ 
ation  is  dependent  on  preexisting  ice  particles,  is  examined. 
The  performance  of  the  microphysical  model  in  the  different 
microphysical  scenarios  and  a  number  of  temperature  sce¬ 
narios  with  and  without  the  influence  of  mountain  waves  is 
tested  through  comparisons  with  lidar  measurements  of  PSCs 
made  from  the  NASA  DC-8  on  23  and  25  January  during  the 
SOLVE/THESEO  2000  campaign  in  the  1999-2000  winter 
and  the  effect  of  mountain  waves  on  local  PSC  production  is 
evaluated  in  the  different  microphysical  scenarios.  Mountain 
waves  are  seen  to  have  a  pronounced  effect  on  the  amount  of 
ice  particles  formed  in  the  simulations.  Quantitative  compar¬ 
isons  of  the  amount  of  solids  seen  in  the  observations  and  the 
amount  of  solids  produced  in  the  simulations  show  the  best 
correspondence  when  NAT  formation  is  allowed  to  take  place 
at  temperatures  above  Tjce.  Mountain  wave-induced  temper¬ 
ature  fluctuations  are  introduced  in  vortex-covering  model 
runs,  extending  the  full  1999-2000  winter  season,  and  the 
effect  of  mountain  waves  on  large-scale  PSC  production  is 
estimated  in  the  different  microphysical  scenarios.  It  is  seen 
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that  regardless  of  the  choice  of  microphysics  ice  particles 
only  form  as  a  consequence  of  mountain  waves  whereas  NAT 
particles  form  readily  as  a  consequence  of  the  synoptic  con¬ 
ditions  alone  if  NAT  nucleation  above  T,ce  is  included  in  the 
simulations.  Regardless  of  the  choice  of  microphysics,  the 
inclusion  of  mountain  waves  increases  the  amount  of  NAT 
particles  by  as  much  as  10%.  For  a  given  temperature  sce¬ 
nario  the  choice  of  NAT  nucleation  mechanism  may  alter  the 
amount  of  NAT  substantially;  three-fold  increases  are  easily 
found  when  switching  from  the  scenario  which  requires  pre¬ 
existing  ice  particles  in  order  for  NAT  to  form  to  the  scenario 
where  NAT  forms  independently  of  ice. 


1  Introduction 

Polar  stratospheric  clouds  (PSCs)  are  known  to  be  of  vi¬ 
tal  importance  to  ozone-depleting  processes  in  the  polar  re¬ 
gions  (Tolbert  and  Toon,  2001).  In  particular,  the  formation 
of  solid  PSC  particles  (ice  and  nitric  acid  trihydrate  (NAT) 
particles)  is  important  since  these  particles  may  grow  large 
enough  for  gravitational  sedimentation  to  have  an  impact, 
thereby  depleting  the  stratosphere  of  water  and  nitrogen  com¬ 
pounds.  The  removal  of  these  compounds  serves  to  enhance 
the  ozone-depleting  conditions  (WMO,  1999). 

In  the  northern  polar  stratosphere  synoptic  temperatures 
usually  hover  around  the  existence  temperature  for  NAT 
(WMO,  1999).  Consequently,  the  additional  cooling  caused 
by  mountain  waves  may  be  of  importance  to  the  produc¬ 
tion  of  solid  PSC  types  in  the  northern  polar  vortex.  Moun¬ 
tain  waves  are  known  to  cause  the  formation  of  solid  PSC 
types  (Carslaw  et  al.,  1998a,b)  and  the  locally  produced  solid 
PSCs  may  exist  at  vast  distances  downstream  (Carslaw  et  al., 
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1999).  The  presence  of  NAT  particles  in  the  Arctic  strato¬ 
sphere  does  not  appear  to  be  sporadic;  in  a  multi-year  study 
of  lidar  data  from  Ny  Alesund,  Biele  et  al.  (2001)  found 
wide-spread  presence  of  NAT  particles  interspersed  with  liq¬ 
uid  aerosol.  In  order  to  be  able  to  produce  reasonable  es¬ 
timates  of  the  load  of  solid  type  PSCs  in  the  Arctic,  it  is 
therefore  important  to  be  able  to  evaluate  the  amount  of  PSCs 
produced  in  mountain  waves  compared  to  the  amount  of  syn- 
optically  produced  PSCs.  Lidar  observations  are  an  excellent 
tool  for  studying  local-scale  PSC  phenomena  in  great  detail 
(Toon  et  al.,  2000).  Detailed  studies  of  lidar  observations  of 
mountain  wave  PSCs  combined  with  microphysical  and  op¬ 
tical  modelling  have  been  carried  out  by  e.g.  Hu  et  al.  (2002) 
and  Fueglistaler  et  al.  (2003). 

According  to  a  study  of  mountain  wave  PSCs  by 
Fueglistaler  et  al.  (2003),  ECMWF  trajectories  lack  the  abil¬ 
ity  to  reproduce  the  cooling  rates  observed  in  mountain 
waves.  In  the  present  study,  a  method  for  incorporating 
mountain  wave  effects  in  microphysical  PSC  simulations  on 
local  as  well  as  synoptic  scales  has  been  developed,  thereby 
allowing  for  an  estimate  of  the  importance  of  mountain  wave 
effects  on  the  PSC  load. 

Whereas  the  existence  of  NAT  particles  in  polar  strato¬ 
spheric  clouds  is  well  established  (Voigt  et  al.,  2000),  the 
formation  mechanism  of  NAT  is  still  debated.  NAT  is  be¬ 
lieved  to  form  by  heterogeneous  nucleation  on  pre-existing 
ice  particles,  requiring  temperatures  below  Tjce  at  some  point 
(Carslaw  et  al.,  1998a,b),  but  recent  studies  have  indicated 
the  need  for  NAT  nucleation  mechanisms  which  are  indepen¬ 
dent  of  the  existence  of  ice  particles,  only  requiring  temper¬ 
atures  below  the  existence  temperature  for  NAT,  a  condition 
which  is  easier  to  meet  than  the  demand  for  pre-existing  ice 
particles  (Drdla  et  al.,  2002;  Pagan  et  al.,  2004;  Irie  et  al., 
2004).  Here,  the  effect  of  a  proposed  nucleation  mechanism 
for  NAT  which  is  active  above  Tlce  is  examined  on  both  local 
and  synoptic  scales. 

In  Sect.  2  the  microphysical  model  used  for  this  study  is 
presented  and  the  method  for  including  mountain  wave  ef¬ 
fects  in  the  simulations  is  introduced  in  Sect.  2.1.  Next,  com¬ 
parisons  between  lidar  measurements  of  PSCs  and  model 
runs  with  variable  microphysics  and  with  and  without  moun¬ 
tain  waves  are  presented  in  Sect.  3.  Finally,  in  Sect.  4  the 
influence  of  mountain  waves  and  the  choice  of  microphysics 
on  the  production  of  solid  PSCs  inside  the  vortex  is  evaluated 
for  the  winter  1999-2000. 


2  The  microphysical  model 

Simulations  of  the  PSC  production  have  been  made  using  the 
Danish  Meteorological  Institute  microphysical  PSC  model 
(Larsen,  1991,  2000).  The  core  of  the  PSC  model  is  a  box 
model  which  for  a  given  time  step  calculates  the  changes  in 
size  distributions,  chemical  composition,  and  physical  phase 
of  an  ensemble  of  particle  types.  The  input  to  the  box  model 


at  each  time  step  is  the  temperature,  pressure,  partial  pres¬ 
sures  of  H?0  and  HNO3,  the  current  number  densities,  and 
the  amount  of  bound  acid  and  water  in  each  size  bin  of  the 
different  aerosol  types,  and  the  box  model  returns  new  val¬ 
ues  of  these  parameters  in  each  size  bin.  Mixtures  of  four 
different  particle  types  are  recognised  by  the  model:  1)  Liq¬ 
uid  particles  assumed  to  be  sulphate  aerosols  turning  into  su¬ 
percooled  ternary  H2O-H2SO4-HNO3  solutions  (STS)  at  low 
temperatures  through  uptake  of  HNO3  and  H2O,  2)  frozen 
sulphate  aerosols  assumed  to  be  sulphuric  acid  tetrahydrate, 
SAT,  3)  NAT  particles  with  inclusions  of  SAT,  and  4)  ice  par¬ 
ticles  with  inclusions  of  NAT  and  SAT.  Each  of  the  four  dif¬ 
ferent  particle  types  has  its  own  size  distribution.  The  size 
distributions  are  discretized  into  N  size  bins  on  a  geometric 
volume  scale.  Within  a  given  size  distribution  particles  are 
shifted  to  higher  or  lower  radius  bins  through  vapour  depo¬ 
sition  and  evaporation.  When  phase  transitions  occur  (e.g. 
all  NAT  evaporating  to  form  a  SAT  particle)  particles  are 
transferred  from  one  size  distribution  to  another.  During  this 
transfer  the  particles  are  assumed  to  have  an  unaltered  radius. 
When  calculating  the  evaporation  and  condensational  growth 
of  liquid,  NAT,  and  ice  particles,  vapour  pressures  over  STS 
are  taken  from  Luo  et  al.  (1995),  vapour  pressures  over  NAT 
are  taken  from  Hanson  and  Mauersberger  (1988),  and  vapour 
pressures  over  ice  are  taken  from  Marti  and  Mauersberger 
(1993).  Ice  particles  form  by  homogeneous  nucleation  out 
of  STS  solutions  at  temperatures  3^1  K  below  the  ice  frost 
point,  Tice,  and  a  homogeneous  nucleation  rate,  derived  from 
experimental  data,  is  used  in  the  simulations  (Koop  et  al., 
2000). 

Recently,  studies  have  been  published  indicating  the  need 
for  a  freezing  process  active  above  7;,r  in  order  to  explain 
observations  of  solid  PSC  particles  (Drdla  et  al.,  2002;  Pa¬ 
gan  et  al.,  2004;  Irie  et  al.,  2004;  Larsen  et  al.,  2004).  As 
an  example  of  a  NAT  nucleation  mechanism  capable  of  gen¬ 
erating  NAT  at  temperatures  above  the  ice  frost  point,  NAT 
nucleation  in  the  present  study  is  described  by  the  nucleation 
rate  of  nitric  acid  dihydrate,  NAD,  found  in  Tabazadeh  et  al. 
(2002)  (assuming  an  instantaneous  conversion  of  NAD  par¬ 
ticles  to  NAT),  although  corrected  by  a  factor  of  0. 1  in  order 
to  comply  with  restrictions  posed  by  observations  of  particle 
size  distributions  and  gas  phase  mixing  ratios  of  HNO3  ac¬ 
cording  to  the  findings  in  Larsen  et  al.  (2004).  The  proposed 
nucleation  mechanism  is  currently  debated,  see  e.g.  Knopf 
et  al.  (2002),  but  it  is,  to  our  knowledge,  the  latest  published 
parameterisation  of  NAT  nucleating  out  of  STS  at  tempera¬ 
tures  above  7)re  and  may  as  such  be  seen  as  a  useful  repre¬ 
sentative  of  scenarios  allowing  for  NAT  nucleation  at  temper¬ 
atures  above  Tjce.  Whether  this  particular  nucleation  mecha¬ 
nism  alone  is  responsible  for  any  NAT  formation  at  temper¬ 
atures  above  Tjce  or  whether  other  nucleation  mechanisms, 
e.g.  heterogeneous  nucleation  (Drdla  et  al.,  2002),  are  active 
cannot  be  concluded  from  the  present  study,  which  only  ad¬ 
dresses  the  effect  of  a  single  proposed  nucleation  mechanism 
for  NAT  which  is  active  above  7/Cf>  as  compared  to  a  scenario 
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where  NAT  formation  requires  the  presence  of  ice  particles. 
Note  that  the  choice  of  nucleation  mechanism  for  NAT  will 
affect  the  formation  of  ice  particles.  If  NAT  can  form  at  tem¬ 
peratures  above  the  ice  frost  point,  ice  may  nucleate  hetero¬ 
geneously  on  preexisting  NAT  particles  as  soon  as  temper¬ 
atures  drop  slightly  below  Tjce  instead  of  only  forming  by 
homogeneous  nucleation  at  temperatures  3-4  K  below  Tjce. 

The  model  is  initialised  by  profiles  of  HNO3  and  H?0 
and  a  background  population  of  sulphuric  aerosols.  The  size 
distribution  of  the  sulphuric  aerosols  is  based  on  a  SAGE  I 
and  SAGE  II  (Stratospheric  Aerosol  and  Gas  Experiment) 
surface  area  climatology  (Hitchman  et  al.,  1994).  Aerosol 
extinctions  are  retrieved  from  the  SAGE/SAM  stratospheric 
aerosol  climatology  Hitchman  et  al.  (1994)  by  interpolation 
in  latitude  and  altitude  in  the  tabulation  of  decadal  mean  val¬ 
ues  and  converted  to  surface  area  densities  by  the  provided 
formula  (Hitchman  et  al.,  1994).  Assuming  the  aerosol  is 
characterised  by  a  lognormal  size  distribution  with  median 
radius  0.0725 /xm  and  geometric  standard  deviation  1.86, 
Pinnick  et  al.  (1976),  the  surface  area  density  is  converted  to 
aerosol  number  density,  resulting  in  values  in  the  order  of  10 
particles  per  cm3,  consistent  with  measurements  (Schreiner 
et  al.,  2003).  The  HNO3  content  is  initialised  by  a  LIMS 
profile  (Gille  and  Russell,  1984)  and  the  HiO  content  is  de¬ 
scribed  as  a  function  of  the  potential  temperature,  ©,  based 
on  a  series  of  frost  point  hygrometer  measurements  in  the 
northern  polar  vortex  (Ovarlez  et  al.,  2004).  Comparing  the 
LIMS  profile  to  observations  of  the  HNO3  profile  made  dur¬ 
ing  SOLVE/THESEO-2000,  Coffey  et  al.  (2002),  the  LIMS 
profile  is  lower  than  the  observations  from  December  and 
higher  than  the  observations  from  March.  It  is  our  opinion 
that  a  profile  placed  in  between  the  two  winter  extremes  of 
December  and  March  is  a  reasonable  value  to  use  for  the  ini¬ 
tialisation  of  the  model. 

The  input  to  the  model  is  provided  by  temperature  trajec¬ 
tories  based  on  ECMWF  analyses.  The  calculation  of  the 
trajectories  uses  6  hourly  ECMWF  analyses  on  a  1.5°  x  1.5° 
longitude-latitude  grid.  The  trajectories  themselves  are  cal¬ 
culated  on  an  equal-area  grid  with  a  grid  distance  of  139  km 
and  start  inside  or  at  the  edge  of  the  polar  vortex  at  9  standard 
isentropic  levels  (360  K,  380  K,  400  K,  435  K,  475  K,  515  K, 
550  K,  600  K,  and  650  K)  in  the  case  of  the  hemispheric  sim¬ 
ulations  and  on  1 1  levels  (400  K,  425  K,  450  K,  475  K,  500  K, 
530  K,  560  K,  590  K,  600  K,  625  K,  and  700  K)  in  the  case  of 
the  small-scale  simulations.  Diabatic  cooling  is  taken  into 
account  in  the  calculations  (Morcrette,  1991;  Knudsen  and 
Grooss,  2000;  Larsen  et  al.,  2002).  At  each  point  of  a  trajec¬ 
tory  the  temperature,  pressure,  potential  temperature,  poten¬ 
tial  vorticity,  latitude,  and  longitude  are  given.  Depending  on 
the  number  of  trajectories  and  their  duration,  microphysical 
simulations  on  a  variety  of  scales  may  be  performed;  from 
very  detailed  model  runs  closely  matching  the  time  and  lo¬ 
cation  of  different  sets  of  measurements  to  vortex-covering 
simulations  extending  over  entire  winter  seasons. 


2.1  Including  mountain  wave  effects  in  the  simulations 

Temperature  fluctuations  due  to  mountain  waves  may  in¬ 
fluence  the  formation  of  solid  type  PSCs  (Carslaw  et  al., 
1998a).  Climatological  studies  of  the  stratosphere  over  Scan¬ 
dinavia  have  shown  that  practically  all  ice  particle  events 
in  this  area  are  attributable  to  the  presence  of  mountain 
waves  (Dornbrack  and  Leutbecher,  2001).  In  order  to  ob¬ 
tain  reliable  estimates  of  the  PSC  load,  the  inclusion  of 
mountain  wave  effects  in  the  modelling  of  PSCs  is  cru¬ 
cial.  ECMWF  reanalyses  lack  sufficient  resolution  to  resolve 
the  full  spectrum  of  mesoscale  mountain  waves  that  occur 
in  the  stratosphere,  see  e.g.  (Fueglistaler  et  al.,  2003).  To 
provide  estimates  of  mountain  wave  influences  on  the  PSC 
load,  mesoscale  temperature  variances  from  hindcast  runs 
using  Version  2  of  the  Naval  Research  Laboratory  Moun¬ 
tain  Wave  Forecast  Model  (MWFM-2)  through  global  anal¬ 
ysis  fields  issued  by  NASA’s  Global  Modeling  and  Analy¬ 
sis  Office  (GMAO)  during  the  Arctic  winter  of  1999-2000 
are  used  to  supplement  the  ECMWF  temperature  trajecto¬ 
ries.  Details  of  the  MWFM-2  algorithms  are  given  by  Eck- 
ermann  and  Preusse  (1999),  Hertzog  et  al.  (2002)  and  Jiang 
et  al.  (2004).  These  fields  are  issued  as  averaged  gridded 
wave-induced  temperature  variance  fields  on  pressure  sur¬ 
faces,  derived  from  the  raw  mountain  wave  ray  data  gen¬ 
erated  by  the  hindcasts,  to  make  them  more  amenable  for 
use  in  offline  transport  simulations.  The  formation  and  use 
of  this  1999-2000  gridded  product  is  described  in  detail  by 
Pierce  et  al.  (2003)  and  Pagan  et  al.  (2004).  An  important 
point  to  note  is  that  the  procedure  of  avaraging  into  a  grid¬ 
ded  variance  product  tends  to  significantly  suppress  typical 
intra-gridbox  variability  found  within  unaveraged  mountain 
wave  fields  (Hertzog  et  al.,  2002;  Brogniez  et  al.,  2003)  and 
thus  the  data  product  used  in  this  study  should  be  viewed 
as  a  working  lower  bound  on  likely  mesoscale  temperature 
variability  due  to  mountain  waves.  Thus,  any  microphysical 
changes  caused  by  insertion  of  these  mountain  wave  fields 
are  likely  significant. 

In  this  study,  the  MWFM-2  fields  are  used  to  introduce 
mountain  wave-induced  temperature  fluctuations  into  the 
simulations.  The  MWFM-2  hindcasts  and  averaging  produce 
daily  fields  of  mountain  wave -induced  temperature  fluctua¬ 
tions  on  a  l°xl°  grid,  specifying  the  amplitude  of  the  tem¬ 
perature  fluctuation  at  seven  standard  pressure  levels  (10,  20, 
30,  40,  50,  70,  and  100  hPa).  Through  interpolations  in  pres¬ 
sure  and  time  the  value  of  the  temperature  amplitude  is  found 
at  the  nearest  grid  point  of  each  trajectory  position,  thereby 
producing  a  data  base  of  mountain  wave-induced  tempera¬ 
ture  fluctuations  at  the  times  and  locations  of  the  trajectories. 
When  running  the  microphysical  simulations,  different  tem¬ 
perature  scenarios  with  and  without  the  inclusion  of  moun¬ 
tain  wave  effects  can  be  examined:  Case  1:  T—Tq,  Case 
2:  T—To—Ta,  Case  3:  T—Tq— T\  —  Tcorr,  where  T  is  the 
temperature  used  as  input  in  the  microphysical  simulations. 
To  is  the  trajectory  temperature  according  to  the  ECMWF 
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Fig.  1.  Flight  tracks  of  the  DC-8  flights  of  23  and  25  January  2000, 
respectively. 


Table  1.  Threshold  values  for  type  classification  of  DC-8  lidar  mea¬ 
surements.  B:  Backscatter  ratio,  D:  Depolarisation. 


DC-8  lidar  threshold  values 

B 

D(%) 

PSC  lb 

0.18  <B 

D<  2.5 

PSC  la 

0.18<B<5.0 

D>  2.5 

PSC  2 

B>  5.0 

D>  2.5 

analysis,  Ta  is  the  mountain  wave-induced  temperature  am¬ 
plitude,  and  TCorr  is  an  additional  temperature  correction, 
based  on  temperature  restrictions  posed  by  the  presence  of 
ice  particles  in  observations.  This  additional  correction  will 
be  described  in  further  detail  below.  Subtraction  of  the  am¬ 
plitude  of  the  mountain  wave  perturbation  in  cases  2  and  3 
corresponds  to  the  maximal  lowering  of  the  temperatures  and 
one  should  therefore  interpret  these  results  as  lower  limits  of 
the  temperature. 


3  Comparing  model  runs  and  lidar  measurements 

Given  the  high  altitude  and  time  resolution  of  lidar  measure¬ 
ments,  very  detailed  information  about  the  observed  aerosols 
can  be  provided  by  such  measurements.  In  order  to  make 
comparisons  between  observations  and  model  results,  it  is 
important  to  ensure  that  the  model  results  closely  match  the 
time  and  location  of  the  measurements.  In  the  present  study, 
lidar  data  from  the  NASA  DC-8  flights  of  23  and  25  Jan¬ 
uary  2000,  during  the  SOLVE/THESEO  2000  campaign  have 
been  considered.  Data  from  these  flights  have  been  studied 
extensively  by  Hu  et  al.  (2002)  and  Fueglistaler  et  al.  (2003). 
A  summary  of  the  SOLVE/THESEO  2000  campaign  can  be 
found  in  Newman  et  al.  (2002).  A  description  of  the  lidar 
system  used  for  the  measurements  can  be  found  in  Browell 
(1989)  and  Browell  et  al.  (1990).  A  map  showing  the  flight 
tracks  of  the  two  flights  is  shown  in  Fig.  1 .  Sets  of  back  tra¬ 
jectories  have  been  initiated  along  the  flight  tracks.  Running 
these  trajectories  produces  information  about  the  aerosol  at 
the  position  and  time  of  the  lidar  measurements  according  to 
the  microphysical  model. 

The  accuracy  of  the  temperatures  and  the  mountain  wave 
corrections  can  be  evaluated  by  considering  the  restrictions 
on  the  temperatures  posed  by  the  presence  of  ice  particles. 
In  areas  where  ice  is  observed,  the  temperature  must  nec¬ 
essarily  be  equal  to  or  colder  than  the  ice  frost  point,  Tice. 
By  keeping  track  of  the  times  and  altitudes  where  ice  parti¬ 
cles  are  observed  one  may  compare  the  model  temperatures 
in  these  areas  to  Tjce.  In  the  DC-8  data,  aerosol  backscatter 
ratios  at  603  nm  greater  than  5.0  with  enhanced  depolarisa¬ 
tion  (>2.5%)  are  considered  to  be  ice  particles,  see  Table  1. 
Comparisons  of  the  model  temperatures  and  Tjce  indicate  a 
warm-bias  in  the  model  temperatures  and  an  additional  tem¬ 
perature  correction,  Tcorr,  is  thus  introduced.  Tcorr  is  defined 
as  the  average  difference  between  the  model  temperature  and 
Tice  in  those  areas  where  ice  particles  are  observed.  It  is  as¬ 
sumed  that  the  additional  temperature  correction  is  associ¬ 
ated  with  an  under-estimation  of  the  mountain  wave-induced 
temperature  fluctuations  and  Tcorr  is  only  introduced  in  the 
simulations  if  the  mountain  wave-induced  temperature  am¬ 
plitude  Ta  exceeds  a  value  of  0.5  K.  The  value  of  Tcorr  is 
found  to  be  1.16K  for  23  January  2000,  and  3.70K  for  25 
January  2000. 

The  model  allows  the  calculation  of  the  backscatter  ratio 
and  depolarisation  based  on  a  T-matrix  calculation  code  as 
described  in  Larsen  et  al.  (2004).  A  number  of  assumptions 
about  the  shape  of  the  particles  and  their  refractive  indices 
must  be  made  and  the  calculated  backscatter  ratio  and  depo¬ 
larisation  will  be  rather  sensitive  to  the  choice  of  these  pa¬ 
rameters.  Hence,  the  calculated  optical  variables  are  only 
suited  for  general  qualitative  comparisons  with  the  measured 
data  and  not  for  detailed,  quantitative  comparisons.  The  dif¬ 
ference  in  temporal  and  altitudinal  resolution  between  the 
measured  and  modelled  quantities  as  well  as  the  limited  res¬ 
olution  (l°xl°)  of  the  mountain  wave-induced  temperature 
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Fig.  2.  Observed  and  modelled  PSC  types  for  23  January,  2000.  Grey  areas  in  the  plots  of  measured  data  represent  valid  data  with  backscatter 
values  less  than  the  minimum  threshold  value  for  PSC  presence  whereas  white  areas  indicate  lack  of  valid  data.  Panel  (A)  shows  the  observed 
PSC  types,  panel  (B)  shows  the  PSC  types  according  to  the  microphysical  model  in  the  T=Tq  scenario,  panel  (C)  shows  the  modelled  PSC 
types  in  the  T=Tq—T^  scenario,  and  panel  (D)  shows  the  modelled  PSC  types  in  the  T =Tf)—T^—Tcorr  scenario. 


corrections  point  to  overall  qualitative  comparisons  instead 
of  detailed  point-by-point  comparisons.  Therefore,  com¬ 
parisons  between  observational  data  and  modelled  data  are 
based  on  cloud  types  instead  of  a  direct  comparison  of  the 
optical  variables.  In  the  case  of  the  observations,  type  dis¬ 
tinction  is  made  by  considering  the  backscatter  ratio  and  the 
depolarisation.  Such  a  cloud  type  distinction  is  only  possible 
when  both  the  backscatter  ratio  and  the  depolarisation  are 
available  and,  hence,  the  type  plots  presented  here  only  rep¬ 
resent  the  areas  where  both  types  of  data  are  present.  Points 
where  the  backscatter  ratio  B  lies  in  the  interval  0.18<B<5.0 
and  the  depolarisation  D  is  larger  than  2.5%  are  classified  as 
type  la  PSCs  whereas  points  where  5.0 <B  and  2.5%  <  D  are 
classified  as  type  2  PSCs  and  points  where  the  depolarisa¬ 
tion  is  smaller  than  2.5%  and  the  backscatter  ratio  is  higher 
than  0.18  (0.18<,B  and  D< 2.5%)  are  considered  to  be  type 
lb  PSCs.  An  overview  of  the  PSC  type  classification  scheme 
for  the  observational  data  is  given  in  Table  1.  It  should  be 
noted  that  in  the  present  paper  only  the  aerosol  backscatter 
ratio  and  not  the  total  backscatter  ratio  is  considered. 


Observations  of  PSCs  indicate  that  very  often  there  are 
both  liquid  and  solid  particles  present  in  the  cloud  (Biele 
et  al.,  2001;  Larsen  et  al.,  2004).  Whether  a  given  PSC  obser¬ 
vation  will  be  classified  as  a  type  la  or  type  lb  PSC  depends 
on  whether  NAT  particles  or  STS  particles  dominate  the  opti¬ 
cal  signal.  When  determining  the  cloud  type  of  a  given  mod¬ 
elled  particle  distribution  it  is  therefore  necessary  to  have  an 
estimate  of  the  proportions  in  the  mixture  of  liquid  and  solid 
particles  which  gives  rise  to  a  classification  as  a  type  la  or 
type  lb  PSC.  In  simulations  of  the  volume  density  of  STS 
and  NAT  particles  and  the  total  colour  index  and  extinction 
a  clear  switch  in  the  colour  index  and  extinction  from  what 
would  be  considered  a  type  lb  cloud  to  a  type  la  cloud  took 
place  when  the  surface  area  density  of  NAT  was  approxi¬ 
mately  5%  of  the  surface  area  density  of  STS,  see  Fig.  5  in 
(Larsen  et  al.,  2004).  When  classifying  the  PSC  type  of  a 
given  model  data  point  this  is  used  as  a  threshold  value  to 
separate  PSC  la  from  PSC  lb;  if  the  surface  area  density  of 
either  NAT  or  ice  in  a  given  point  constitutes  5%  or  more  of 
the  surface  area  density  of  STS  particles,  the  point  is  classi¬ 
fied  as  a  type  1  a  or  type  2  point,  respectively.  If  not,  the  point 
is  classified  as  a  type  lb  point. 
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Fig.  3.  Observed  and  modelled  PSC  types  for  25  January  2000.  Grey  areas  in  the  plots  of  measured  data  represent  valid  data  with  backscatter 
values  less  than  the  minimum  threshold  value  for  PSC  presence  whereas  white  areas  indicate  lack  of  valid  data.  Panel  (A)  shows  the  observed 
PSC  types,  panel  (B)  shows  the  PSC  types  according  to  the  microphysical  model  in  the  T=Tq  scenario,  panel  (C)  shows  the  modelled  PSC 
types  in  the  T =Tq—T2 4  scenario,  and  panel  (D)  shows  the  modelled  PSC  types  in  the  T =TQ—T^  —  Tcorr  scenario. 


Comparisons  of  the  PSC  type  distribution  in  the  measured 
and  the  modelled  data  are  shown  in  Fig.  2  for  23  January 
2000,  and  Fig.  3  for  25  January  2000.  NAT  nucleation  above 
Tjce  has  been  included  in  the  simulations.  In  the  figures,  the 
upper  left  panel  (A)  shows  the  PSC  type  distribution  in  the 
measurement  data.  The  upper  right  panel  (B)  shows  model 
results  from  the  T—Tq  scenario,  the  lower  right  panel  (C) 
shows  the  results  from  the  T—Tq—Ta  scenario  and  the  re¬ 
sults  from  the  T=Tq—Ta~ Tcorr  scenario  are  shown  in  the 
lower  left  panel  (D).  In  the  simulations  of  the  two  DC-8 
flights  NAT  nucleation  above  Tjce  has  been  included. 

In  the  case  of  23  January  2000  (see  Fig.  2),  areas  of  all 
three  PSC  types  are  seen  in  the  observations.  When  con¬ 
sidering  the  model  data  it  is  seen  that  in  the  case  where  no 
temperature  corrections  have  been  introduced  (panel  B)  the 
model  handles  the  presence  of  PSC  la  and  PSC  lb  reason¬ 
ably  well.  Around  09:00  UT  the  observations  indicate  PSC 
la  around  21-23  km  followed  by  PSC  lb  from  19  to  21  km. 
In  the  model,  this  region  shows  both  PSC  la  and  PSC  lb  as 
well,  although  the  PSC  la  field  extends  too  far  downwards 
and  the  PSC  lb  field  is  lying  too  low.  Around  10:00  UT 


the  observations  indicate  an  area  of  no  PSCs  between  15  and 
19  km  and  from  22  to  23  km,  possibly  with  small  amounts 
of  PSC  la  clouds  in  the  lower  area.  No  valid  data  is  present 
covering  the  intermediate  altitude  interval.  In  the  model,  the 
entire  altitude  interval  contains  PSCs,  with  a  small  area  of 
PSC  la  at  the  top  of  the  profile  and  PSC  lb  in  the  interme¬ 
diate  and  lower  part.  The  next  section  in  the  observational 
data  stretches  from  12:00  to  13:30  UT  in  the  altitude  inter¬ 
val  from  16  to  24  km  and  shows  no  PSCs  in  the  early  and/or 
lower  parts  of  the  segment  and  PSC  la  in  the  higher  part.  In 
the  model  data  PSC  la  is  the  prominent  feature,  although  the 
PSCs  are  seen  to  extend  too  low  and  span  a  longer  time  inter¬ 
val  than  the  observations.  The  last  segment  of  the  observa¬ 
tions  show  distinct  areas  of  PSC  2  (from  22  to  24  km  between 
15:30  and  16:00  UT),  PSC  la  (from  15:00  to  17:00  UT  and 
from  19  km  and  upwards)  and  PSC  lb  which  is  mainly  seen 
from  16:00  to  17:00  UT  between  17  and  23  km.  In  the  model 
data,  the  PSC  2  field  is  not  reproduced.  Ample  amounts  of 
PSC  lb  are  observed  as  well  as  PSC  la,  although  the  PSC 
la  field  is  generally  too  small.  Overall,  the  model  data  ap¬ 
pears  to  extend  the  PSC  fields  too  far  downwards,  which  may 
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Fig.  4.  Comparison  of  measured  and  modelled  PSC  types  for  23 
January  2000.  Panel  (A)  shows  the  observed  PSC  types  whereas 
panel  (B)  shows  the  modelled  PSC  types  in  the  microphysical  sce¬ 
nario  where  NAT  nucleation  above  the  ice  frost  point  is  included 
and  panel  (C)  shows  the  modelled  PSC  types  in  the  microphysical 
scenario  where  NAT  nucleation  above  the  ice  frost  point  is  not  in¬ 
cluded.  In  both  of  the  model  simulations  the  coldest  temperature 
scenario  T =T$—T &  —  Tcorr  has  been  used. 


Fig.  5.  Comparison  of  measured  and  modelled  PSC  types  for  25 
January  2000.  Panel  (A)  shows  the  observed  PSC  types  whereas 
panel  (B)  shows  the  modelled  PSC  types  in  the  microphysical  sce¬ 
nario  where  NAT  nucleation  above  the  ice  frost  point  is  included 
and  panel  (C)  shows  the  modelled  PSC  types  in  the  microphysical 
scenario  where  NAT  nucleation  above  the  ice  frost  point  is  not  in¬ 
cluded.  In  both  of  the  model  simulations  the  coldest  temperature 
scenario  T =  7b— 7^  — 7corr  has  been  used. 


be  indicative  of  an  altitude  bias  in  the  model  data.  It  should 
be  mentioned,  however,  that  the  altitude  is  determined  from 
the  pressure  of  the  input  trajectories  and  the  scale  height,  and 
should  be  viewed  as  an  approximation.  Most  notably  is  the 
lack  of  ice  particles  in  the  model  data  considering  the  rather 
prominent  ice  feature  seen  in  the  observations.  When  moun¬ 
tain  wave  effects  are  included  in  the  simulations  (panels  C 
and  D)  an  increase  in  PSC  la  is  seen  at  the  expense  of  PSC 
lb.  Ice  clouds  are  only  seen  in  the  coldest  of  the  three  temper¬ 
ature  scenarios,  (panel  D).  The  PSC  la  field  observed  around 
15:00  UT  is  not  reproduced  by  the  model  in  any  of  the  three 
temperature  scenarios. 

Considering  the  DC-8  flight  of  25  January  2000,  all  three 
PSC  types  are  observed  as  well  (see  Fig.  3,  panel  A).  The  first 
section  of  the  observational  data  (from  08:00  to  10:00  UT) 
shows  only  PSC  la  and  PSC  lb,  with  PSC  la  dominating  the 
early  part  of  this  section  and  PSC  lb  dominating  the  late  part. 
In  the  case  of  the  model  data  the  synoptic  scenario  (panel  B) 
shows  a  mixture  of  PSC  la  and  PSC  lb  as  well,  although 
the  PSC  la  field  extends  a  little  later  than  what  is  seen  in 
the  observations.  The  next  section  of  the  observational  data 
from  12:30  to  17:00  UT  shows  three  prominent  PSC  2  areas 


surrounded  by  PSC  la.  The  middle  part  of  this  section  (from 
14:00  to  15:30  UT)  is  a  complex  mixture  of  PSC  la  and  PSC 
lb.  In  the  synoptic  model  data  (panel  B)  the  PSC  la  fields  are 
reproduced  reasonably  well.  Most  notable  is  the  complete 
lack  of  PSC  2  fields  in  the  synoptic  model  data.  Proceeding 
to  scenarios  where  mountain  wave  effects  are  included  (panel 
C  and  D)  an  increase  in  the  amount  of  PSC  la  is  seen  at 
the  expense  of  PSC  lb.  Going  to  the  coldest  of  the  three 
temperature  scenarios,  T—Tq— Ta~ Tcorr  (panel  D)  areas  of 
PSC  2  are  seen  in  the  model  data  as  well.  However,  the  PSC 
2  areas  in  the  latest  segment  span  time  intervals  which  are 
slightly  long  compared  to  the  observations  and  in  the  first 
segment  (around  08:00  UT)  a  PSC  2  field  is  seen  in  the  model 
data  which  is  not  seen  in  the  observations.  Given  that  the 
amount  of  PSC  2  in  the  T —Tq—T A  —  Tcorr  seems  to  be  rather 
large  compared  to  the  observed  amount  of  PSC  2,  one  may 
suspect  that  the  value  of  Tcorr  has  been  too  high  in  this  case. 
Tcorr  is  an  average  value  based  on  the  difference  between 
Tice  and  T  in  the  areas  where  PSC  2  observations  are  made 
and  should  only  be  considered  a  rough  estimate. 
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Fig.  6.  Percentage  of  PSC  la  in  observational  and  model  data  for 
23  January  2000,  for  three  different  temperature  scenarios.  Model 
results  are  shown  for  two  different  microphysical  scenarios,  one 
where  NAT  nucleation  above  Tjce  is  allowed,  and  one  where  this 
is  not  the  case.  Horizontal  black  lines  mark  the  PSC  la  percent¬ 
age  in  simulations  where  the  HNO3  content  has  been  changed  by 
ilppbv.  The  variation  in  the  case  where  NAT  nucleation  above 
Tice  is  not  included  is  not  discernable  on  this  scale. 
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Fig.  7.  Percentage  of  PSC  2  in  observational  and  model  data  for 
23  January  2000,  for  three  different  temperature  scenarios.  Model 
results  are  shown  for  two  different  microphysical  scenarios,  one 
where  NAT  nucleation  above  Tjce  is  allowed,  and  one  where  this 
is  not  the  case.  The  variation  in  the  percentage  of  PSC  2  when  vary¬ 
ing  the  HNO3  content  by  ±  1  ppbv  is  not  discernable  on  this  scale. 


Overall,  it  seems  that  the  microphysical  model  is  capable 
of  producing  PSC  type  distributions  which  correspond  rea¬ 
sonably  well  with  observations.  In  both  the  examined  cases 


the  synoptic  temperature  scenario  ( T —Tq)  was  unable  to  re¬ 
solve  the  observed  type  2  PSCs.  Going  to  mountain  wave- 
influenced  scenarios  resulted  in  PSC  2  presence  at  the  correct 
locations,  even  though  indications  of  an  overestimated  value 
of  Tcorr  was  seen  in  the  case  of  25  January  2000.  As  a  whole, 
the  best  agreement  between  the  observations  and  the  corre¬ 
sponding  simulations  is  seen  in  the  cases  where  the  influence 
of  mountain  waves  has  been  included  since  these  cases  seem 
to  result  in  a  better  reproduction  of  the  solid  PSC  types. 

The  simulations  used  for  comparison  with  the  observa¬ 
tion  data  above  all  included  NAT  nucleation  at  temperatures 
above  7}ce.  In  order  to  illustrate  the  difference  between  the 
outcome  of  simulations  with  and  without  this  particular  nu¬ 
cleation  mechanism  comparisons  between  the  two  different 
microphysical  scenarios  are  shown  in  Fig.  4  for  23  January 
2000,  and  in  Fig.  5  for  25  January  2000.  In  both  of  these  fig¬ 
ures,  panel  (A)  shows  the  observed  distribution  of  PSC  types 
whereas  panel  (B)  shows  the  modelled  PSC  distribution 
when  NAT  nucleation  above  Tjce  is  included  in  the  simula¬ 
tions  and  panel  C  shows  the  modelled  PSC  distribution  when 
this  particular  nucleation  mechanism  is  not  included  in  the 
simulations  and  NAT  forms  on  pre-existing  ice  particles  in¬ 
stead.  The  coldest  temperature  scenario,  T =Tq—T \—TCOrr, 
is  used  in  all  the  simulations. 

In  the  case  of  23  January  2000  (Fig.  4),  large  differences 
are  seen  between  the  two  different  microphysical  scenarios; 
when  NAT  nucleation  above  Tjce  is  allowed  (panel  B)  large 
sections  of  PSC  la  are  seen,  largely  corresponding  to  the 
PSC  la  fields  seen  in  the  observations  (panel  A).  A  stretch 
of  PSC  2  is  seen  as  well,  although  slightly  delayed  compared 
to  the  observations.  Going  to  the  scenario  where  NAT  is  only 
allowed  to  nucleate  on  pre-existing  ice  particles  the  model  re¬ 
sults  are  completely  dominated  by  PSC  lb  and  only  sporadic 
traces  of  PSC  la  are  seen.  This  is  not  in  correspondence  with 
the  observations  which  show  ample  amounts  of  PSC  la. 

The  comparison  for  25  January  2000,  is  shown  in  Fig.  5. 
The  simulation  where  NAT  only  forms  by  nucleation  on 
pre-existing  ice  particles  (panel  C)  is  dominated  by  PSC  lb 
(panel  C),  but  shows  PSC  2  fields  in  the  correct  positions. 
Compared  to  the  measurements  (panel  A)  the  amount  of  PSC 
la  is  underestimated  in  this  scenario.  In  the  case  where  NAT 
nucleation  above  Tjce  is  allowed  (panel  B)  much  more  pro¬ 
nounced  fields  of  PSC  la  are  seen,  which  corresponds  better 
with  the  observations  than  the  case  where  NAT  nucleation 
above  7}ce  is  not  allowed  (panel  C).  Both  microphysical  sce¬ 
narios  display  a  PSC  2  area  around  08:00  UT  which  is  not 
seen  in  the  observations,  an  indication  that  the  final  tempera¬ 
ture  correction  Tcorr  may  be  over-estimated  in  the  case  of  25 
January  2000. 

In  the  qualitative  comparisons  between  observations  and 
simulations  described  above  indications  of  a  preference  for 
mountain  wave-perturbed  scenarios  where  NAT  nucleation 
above  T/ce  is  allowed  was  found.  In  the  following,  this  qual¬ 
itative  analysis  is  extended  to  a  quantitative  one  based  on 
the  overall  PSC  statistics  of  the  entire  DC-8  flights  and  the 
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corresponding  model  simulations.  For  the  sections  of  the 
model  data  which  correspond  to  areas  where  both  backscat- 
ter  ratio  and  depolarisation  data  are  available  and,  hence,  a 
PSC  type  classification  is  possible,  it  is  investigated  whether 
or  not  any  solid  PSCs  are  present  in  the  model  data  and  the 
percentage  of  model  points  containing  either  PSC  la  or  PSC 
2  can  be  determined.  The  percentage  of  PSC  2  and  PSC  la 
in  the  observational  data  (red  columns)  and  in  the  simula¬ 
tions  is  shown  in  Fig.  6  and  Fig.  7  for  23  January  2000,  and 
in  Fig.  8  and  Fig.  9  for  25  January  2000.  In  the  figures,  the 
percentages  of  PSC  la  and  PSC  2  in  the  model  data  are  deter¬ 
mined  for  two  sets  of  model  data:  one  where  NAT  nucleation 
above  7)ce  is  allowed  as  described  earlier  in  this  paper  (blue 
columns),  and  one  where  NAT  nucleation  above  7/C(?  is  not 
allowed  (green  columns).  In  the  latter  case,  NAT  forms  on 
preexisting  ice  particles. 

For  both  23  and  25  January  2000,  the  percentage  of  PSC  la 
in  the  model  data  shows  the  best  correspondence  with  the 
measurements  when  NAT  nucleation  above  Tjce  is  included 
(see  Fig.  6  and  Fig.  8)  than  when  this  is  not  the  case.  When 
NAT  nucleation  above  7/,.,,  is  included  in  the  microphysics 
the  percentage  of  PSC  la  in  the  model  data  is  of  the  same 
order  of  magnitude  as  the  percentage  of  PSC  la  in  the  obser¬ 
vations  whereas  the  percentage  of  PSC  la  in  the  case  where 
NAT  nucleation  is  not  included  is  much  smaller.  In  the  case 
of  the  PSC  la,  the  influence  of  mountain  waves  is  crucial  in 
the  case  where  NAT  nucleation  above  Tjce  is  not  included  in 
the  microphysics:  In  this  case  PSC  la  is  only  seen  as  a  conse¬ 
quence  of  mountain  waves  whereas  PSC  la  is  seen  in  ample 
amounts  already  as  a  consequence  of  the  synoptic  conditions 
alone  when  NAT  nucleation  above  Ttce  is  allowed.  For  both 
the  days  considered  here,  the  synoptic  temperatures  alone 
were  sufficient  to  generate  the  majority  of  the  observed  PSC 
la  and  in  such  cold  cases  where  the  synoptic  temperatures  al¬ 
ready  are  below  Tnat  any  additional  mountain  wave  cooling 
will  not  cause  any  large  changes.  In  other  cases,  where  the 
synoptic  temperatures  are  warmer,  mountain  wave  cooling 
may  have  a  larger  influence  on  the  amount  of  NAT  particles 
in  the  scenario  where  NAT  nucleation  is  allowed  above  Tjce. 

In  a  recent  study  it  was  established  that  PSC  la  obser¬ 
vations  from  the  early  part  of  the  winter  1999-2000  could 
not  be  explained  by  the  presence  of  mountain  waves  and  that 
some  NAT  nucleation  mechanism  active  above  the  ice  frost 
point  was  necessary  in  order  to  explain  the  observations  (Pa¬ 
gan  et  al.,  2004).  The  present  work  supports  this  conclusion 
since,  even  in  the  presence  of  mountain  waves,  the  amount 
of  PSC  la  produced  when  there  is  no  NAT  nucleation  above 
Tice  does  not  correspond  to  the  amount  of  PSC  la  seen  in 
observations.  When  NAT  nucleation  above  7)ce  is  active  a 
better  correspondence  between  observational  data  and  model 
data  is  seen. 

In  the  case  of  the  PSC  2,  the  influence  of  mountain  waves 
is  much  larger,  see  Fig.  7  and  Fig.  9.  In  this  case,  no  PSC 
2  is  seen  in  the  model  data  when  the  synoptic  temperatures 
alone  are  taken  into  consideration,  regardless  of  the  choice 
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Fig.  8.  Same  as  Fig.  6,  but  for  25  January  2000. 


Percentage  of  PSC  2,  25/1-2000 

25-i 


20- 


15- 


NO 

0s 

10- 


Fig.  9.  Same  as  Fig.  7,  but  for  25  January  2000.  Horizontal  black 
lines  mark  the  PSC  2  percentage  when  the  HNO3  content  is  varied 
by  ±1  ppbv.  The  variation  in  the  case  where  N AT  nucleation  above 
Tice  is  not  included  is  not  discemable  on  this  scale. 
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of  microphysics.  This  is  in  agreement  with  the  findings  of 
Dornbrack  and  Leutbecher  (2001),  where  a  climatological 
study  indicated  an  almost  complete  dependence  on  mountain 
waves  for  the  PSC  2  production  over  Scandinavia.  When 
mountain  wave  effects  are  included  PSC  2  is  seen  in  the 
model  data  in  the  coldest  scenario  (T=To—TA—Tcorr),  on 
23  January  2000,  but  only  when  NAT  nucleation  above  7)ce 
is  allowed.  On  25  January  2000,  PSC  2  is  observed  in  both 
microphysical  scenarios  in  the  coldest  temperature  scenario, 
although  the  percentage  of  PSC  2  is  too  large,  with  the  largest 
discrepancy  in  the  case  where  NAT  nucleation  above  7/,,,,  is 
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allowed.  The  fact  that  the  model  overestimates  the  percent¬ 
age  of  PSC  2  on  25  January  2000,  may  be  an  indication  that 
the  value  of  Tcorr  is  too  large  in  this  case,  as  mentioned  ear¬ 
lier.  The  influence  of  the  choice  of  NAT  nucleation  mech¬ 
anism  on  the  percentage  of  PSC  2  is  not  surprising;  when 
NAT  particles  are  present,  ice  may  nucleate  heterogeneously 
as  soon  as  small  supersaturations  with  respect  to  ice  occur, 
whereas  much  larger  supersaturations  are  needed  in  order  for 
homogeneous  nucleation  of  ice  particles  to  take  place.  Over¬ 
all,  the  correspondence  between  the  percentage  of  PSC  2  in 
the  model  data  and  in  the  observations  seems  to  be  best  when 
NAT  nucleation  above  Tjce  is  taken  into  consideration,  as  was 
the  case  for  the  PSC  la.  However,  in  the  case  of  the  PSC  2, 
the  impact  of  mountain  waves  is  much  larger  compared  to  the 
PSC  la.  Type  2  clouds  seem  to  form  exclusively  as  a  conse¬ 
quence  of  mountain  wave  effects,  regardless  of  the  choice  of 
microphysics,  whereas  type  la  clouds  are  observed  already 
as  a  consequence  of  the  synoptic  conditions. 

The  sensitivity  of  the  model  results  towards  changes  in 
the  amount  of  HNO3  has  been  examined  by  performing 
model  runs  where  the  initial  LIMS  profile  has  been  shifted 
by  ±1  ppbv  HNO3.  The  results  from  these  sensitivity  runs 
are  marked  in  Figs.  6-9  by  horizontal  black  lines  (not  all 
cases  are  marked;  in  some  cases  the  variation  is  not  discern- 
able  on  the  scale  used  in  the  figures).  The  amount  of  model 
data  points  classified  as  PSC  la  varies  between  the  different 
HNO3  scenarios,  but  the  differences  between  the  different 
microphysical  scenarios  and  temperature  scenarios  are  sta¬ 
ble  under  the  variations  in  the  HNO3  profile.  It  is  seen  from 
Figs.  6-9  that  the  solid  PSC  percentage  can  decrease  if  the 
amount  of  HNO3  is  increased.  This  has  to  do  with  the  way 
solid  PSCs  are  defined  in  the  type  discrimination  algorithm, 
where  a  data  point  is  defined  as  a  PSC  1  a  or  PSC  2  point  if  the 
surface  area  density  of  either  NAT  or  ice  constitutes  5%  or 
more  of  the  surface  area  density  of  STS.  Changing  the  mix¬ 
ing  ratio  of  HNO3  will  affect  both  the  surface  area  density  of 
STS  and  the  solid  particle  types.  In  particular,  raising  Tn at 
by  increasing  the  mixing  ratio  of  HNO3  may  cause  changes 
in  the  surface  area  density  of  NAT  because  the  cooling  rate 
at  the  point  where  the  temperature  drops  below  That  greatly 
affects  the  size  distribution  of  the  nucleating  NAT  particles, 
e.g.  forming  a  few  large  NAT  particles  instead  of  a  lot  of 
small  particles. 

4  The  influence  of  mountain  waves  on  large-scale  PSC 
production 

In  Sect.  3  it  was  shown  that  mountain  waves  may  have  a  sig¬ 
nificant  influence  on  the  amount  of  ice  particles  and  that  the 
inclusion  of  mountain  wave  effects  was  necessary  in  order  to 
be  able  to  produce  amounts  of  ice  PSC  events  comparable 
to  those  seen  in  lidar  observations.  In  addition,  it  was  seen 
that  including  a  NAT  nucleation  mechanism  which  is  active 
above  Tjce  in  the  simulations  resulted  in  a  better  correspon¬ 


dence  between  the  amount  of  NAT  and  ice  particles  in  the 
observations  and  the  model  output.  The  apparent  importance 
of  mountain  waves  on  ice  PSC  production  and  of  the  choice 
of  NAT  nucleation  mechanism  on  production  of  NAT  and  ice 
on  local  scales  naturally  leads  to  the  question  of  the  influence 
of  mountain  waves  and  NAT  nucleation  above  Tjce  on  large- 
scale  PSC  production.  In  order  to  address  this  issue  vortex¬ 
covering  model  runs  extending  from  15  November  1999,  to 
15  March  2000,  have  been  made.  A  set  of  19  140  trajectories, 
distributed  over  9  isentropic  levels,  is  initiated  on  11  Jan¬ 
uary  2000,  and  calculated  backwards  and  forwards  in  time, 
thereby  covering  the  entire  winter  season.  The  model  is  run 
for  three  different  temperature  scenarios,  one  unperturbed  by 
mountain  waves  and  two  with  increasingly  stronger  moun¬ 
tain  wave  corrections  as  described  in  Sect.  2.1.  In  the  case 
of  the  hemispheric  simulation,  a  direct  estimate  of  Tcorr  is 
not  possible  and  instead,  an  average  of  the  two  values  found 
for  23  and  25  January  2000,  is  used  (Tcorr= 2.4  in  the  hemi¬ 
spheric  runs).  The  value  of  Tcorr  in  the  hemispheric  case 
should  be  considered  a  rough  estimate  since  it  is  based  on 
just  two  local  cases.  Also,  given  that  indications  have  been 
found  that  the  calue  of  Tcorr  for  25  January  2000,  is  over¬ 
estimated,  the  hemispheric  Tcorr  value  should  be  considered 
an  upper  estimate.  In  addition  to  the  different  temperature 
scenarios,  two  microphysical  scenarios  are  tested  as  well, 
one  where  NAT  nucleation  above  T,ce  is  included  and  one 
where  this  is  not  the  case. 

Initially,  the  influence  of  mountain  waves  on  the  tempera¬ 
tures  is  examined  by  analysing  the  trajectory  data  and  deter¬ 
mining  the  percentage  of  the  trajectories  where  the  temper¬ 
ature  drops  below  either  T^at  or  Tlce  as  a  function  of  day 
number.  This  provides  an  estimate  of  the  amount  of  trajecto¬ 
ries  where  NAT  or  ice  particles  can  possibly  exist.  In  Fig.  10 
these  percentages  are  shown  for  three  different  model  lay¬ 
ers.  An  ample  amount  of  trajectories  is  seen  to  have  temper¬ 
atures  below  Tv j  even  without  the  presence  of  mountain 
waves  (panels  A,  B,  C,  red  curves).  However,  the  percent¬ 
age  of  trajectories  where  T <T^at  increases  as  mountain 
wave  effects  are  taken  into  consideration  (panels  A,  B,  C, 
black  and  blue  curves),  with  the  highest  trajectory  frequen¬ 
cies  adhering  to  the  strongest  possible  temperature  pertur¬ 
bation  (T=T0-TA-Tcorr,  panels  A,  B,  C,  blue  curves).  In 
the  case  of  trajectories  where  T  <Tice  (panels  D,  E,  F)  only  a 
small  amount  of  trajectories  fulfills  this  criterion  when  moun¬ 
tain  wave  effects  are  not  taken  into  consideration.  In  this 
case,  temperatures  below  T,ce  are  only  observed  briefly  just 
before  day  0  and  just  before  day  30  and  40  (panels  D,  E, 
F,  red  curves).  When  mountain  wave-induced  temperature 
fluctuations  are  included  the  number  of  trajectories  where 
T <Tjce  increases  and  temperatures  below  7)ce  are  seen  over 
extended  periods  of  time  during  the  winter  season,  although 
the  numbers  typically  are  small  (<0.5%),  with  the  highest 
values  reaching  around  4%  (panels  D,  E,  F,  black  and  blue 
curves).  Again,  the  highest  trajectory  frequencies  are  ob¬ 
served  in  the  case  where  the  strongest  possible  temperature 


Atmos.  Chem.  Phys.,  5,  739-753,  2005 


www.atmos-chem-phys.org/acp/5/739/ 


S.  H.  Svendsen  et  al.:  Influence  of  mountain  waves  and  NAT  nucleation  on  PSC  formation 


749 


vP 

g 


I 

! 


8 

S' 


h 


1  <T|ce,  475 

K 

D 

,  1 

1 

! 

1 

iii 

ii 

0,0  I  ""  ,  iP»PP,"  , 

-60  -40  -20  0  20  40  60  80 


35 
30 
25 
20 
15 
10 
5 
0 

-60  -40  -20  0  20  40  60  80 
4,0 

3.5 
3,0 

2.5 

2,0 

1.5 
1,0 
0,5 
0,0 

-60  -40  -20  0  20  40  60  80 

Daynumber  (2000) 


T<Tice,  51! 

K 

E 

L 

Zl 

L 

m 

Fig.  10.  Percentage  of  trajectories  where  the  temperature  is  below  either  T^at  (top  row)  or  Tlce  (bottom  row)  for  three  different  model 
layers,  475  K,  515K,  and  550  K.  Results  are  shown  for  three  different  temperature  scenarios  according  to  the  colour  coding,  one  without 
mountain  wave  pertubations  and  two  where  increasingly  stronger  mountain  wave  pertubations  are  included. 


perturbation  is  employed  (panels  D,  E,  F,  blue  curves).  As 
mentioned  earlier,  these  results  can  only  provide  an  estimate 
of  the  possible  existence  of  NAT  and  ice  particles.  In  order 
to  evaluate  the  actual  load  of  NAT  and  ice  particles  one  must 
examine  the  outcome  of  the  full  microphysical  simulation. 

The  percentage  of  trajectories  containing  NAT  or  ice  par¬ 
ticles  within  the  polar  vortex  as  a  function  of  day  number  ac¬ 
cording  to  the  microphysical  simulation  is  shown  in  Fig.  11. 
Results  are  shown  for  three  different  temperature  scenarios 
and  for  two  different  choices  of  microphysics.  In  the  present 
analysis,  a  trajectory  is  said  to  contain  ice  or  NAT  if  the  vol¬ 
ume  density  of  these  particle  types  exceeds  zero.  This  is 
a  very  liberal  definition  of  the  presence  of  ice  or  NAT  and 
should  be  interpreted  as  an  upper  limit.  In  the  case  of  NAT 
(panels  A,  B,  C),  it  is  seen  that  NAT  forms  readily  as  a  con¬ 
sequence  of  the  synoptic  conditions  alone  when  NAT  nucle¬ 
ation  above  7jce  is  allowed  (panels  A,  B,  C,  solid  red  curves). 
When  this  is  not  the  case  no  NAT  particles  are  seen  (panels 
A,  B,  C,  dotted  red  curves).  When  mountain  wave  effects 
are  included  the  period  during  which  NAT  particles  are  seen 
(panels  A,  B,  C,  black  and  blue  curves)  starts  earlier  in  the 
winter  compared  to  the  case  where  mountain  wave  effects 
were  not  included  (panels  A,  B,  C,  red  curves).  In  the  cases 


where  mountain  wave  effects  are  taken  into  account  NAT  is 
readily  seen  throughout  the  winter,  regardless  of  the  choice 
of  microphysics,  although  the  percentage  of  trajectories  con¬ 
taining  NAT  is  consistently  higher  in  the  case  where  NAT 
nucleation  above  is  included. 

The  importance  of  mountain  wave  effects  for  NAT  forma¬ 
tion  seems  much  larger  in  the  hemispheric  simulations  than 
in  the  local  ones.  This  could  be  because  the  synoptic  temper¬ 
atures  alone  are  generally  well  below  T^at  in  the  local  scale 
simulations.  If  the  additional  cooling  caused  by  mountain 
waves  is  not  sufficient  to  lower  the  temperatures  below  7’,,,,,, 
the  addition  of  mountain  wave  effects  will  not  cause  many 
changes.  Doing  local  scale  model  runs  in  relatively  warmer 
periods  (e.g.  the  early  part  of  the  winter)  may  reveal  a  greater 
impact  of  mountain  waves  on  NAT  formation  on  local  scales 
as  well. 

Considering  the  ice  particles  (panels  D,  E,  F)  it  is  evident 
that,  regardless  of  the  choice  of  microphysics,  practically  no 
ice  particles  form  as  a  consequence  of  the  synoptic  condi¬ 
tions  alone  (panels  D,  E,  F,  red  curves),  as  was  the  case  in 
the  local-scale  model  runs  as  well.  When  mountain  wave  ef¬ 
fects  are  included  the  percentage  of  trajectories  containing 
ice  particles  increases  noticeably  (panels  D,  E,  F,  black  and 
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Fig.  11.  Percentage  of  trajectories  containing  NAT  (left  column)  and  ice  (right  column)  particles  within  the  polar  vortex  as  a  function  of 
daynumber  for  the  winter  1999-2000.  Results  are  shown  for  three  different  model  layers,  475  K,  515  K,  and  550  K,  and  for  three  different 
temperature  scenarios  according  to  the  colour  coding,  one  without  mountain  wave  pertubations  and  two  where  increasingly  stronger  mountain 
wave  pertubations  are  included.  The  solid  lines  are  from  a  series  of  simulations  where  NAT  nucleation  above  Tjce  is  allowed  whereas  the 
dotted  lines  are  from  a  simulation  sequence  where  this  is  not  the  case. 


blue  curves),  and  trajectories  containing  ice  particles  are  ob¬ 
served  throughout  the  winter  season.  When  NAT  nucleation 
above  is  included  (solid  curves)  the  percentage  is  higher 
than  in  the  case  where  NAT  nucleation  above  l)ce  is  not  al¬ 
lowed  (dotted  curves).  These  results  are  in  agreement  with 
the  findings  of  Dornbrack  and  Leutbecher  (2001)  who  found 
that  practically  all  ice  particle  formation  over  Scandinavia 
was  a  consequence  of  mountain  wave  activity. 


As  seen  above  the  effect  of  localised  mountain  wave  ef¬ 
fects  may  extend  beyond  their  local  scale.  In  hemispheric 
model  runs  where  mountain  wave  effects  have  been  included 
NAT  particle  production  is  seen  to  increase  and  practically  all 
the  ice  particles  produced  in  the  simulations  can  be  directly 
attributed  to  the  effect  of  mountain  waves. 
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5  Results  and  discussion 

Detailed  model  runs  matching  lidar  measurements  from  two 
specific  flight  days  (23  and  25  January  2000)  of  the  NASA 
DC-8  during  SOLVE/THESEO  2000  have  been  analysed  and 
the  correspondence  between  model  data  and  observational 
data  has  been  examined.  Considering  a  quantitative  analysis 
of  the  percentage  of  data  points  containing  PSC  la  or  PSC 
2  the  correspondence  is  better  when  NAT  nucleation  above 
Tjce  is  allowed  than  when  this  is  not  the  case.  In  the  case 
of  the  local  scale  model  runs  the  amount  of  PSC  la  is  only 
slightly  affected  by  the  inclusion  of  mountain  wave  effects 
in  the  simulations,  whereas  mountain  waves  have  a  notice¬ 
able  impact  on  the  amount  of  PSC  2.  In  the  latter  case,  the 
correspondence  between  measurements  and  model  data  in¬ 
creases  greatly  when  mountain  waves  are  added.  The  limited 
effect  of  mountain  waves  on  NAT  formation  in  the  two  sets 
of  examined  DC-8  data  may  be  due  to  the  fact  that  the  syn¬ 
optic  temperatures  alone  are  rather  low;  substantial  amounts 
of  NAT  are  seen  as  a  consequence  of  the  synoptic  conditions 
alone.  Any  additional  cooling  is  not  likely  to  result  in  signifi¬ 
cant  changes  unless  the  cooling  results  in  temperatures  below 
Tjce.  In  a  series  of  sensitivity  tests  the  qualitative  differences 
between  the  different  microphysical  scenarios  were  not  af¬ 
fected  by  variations  in  the  initial  HNO3  profile  of  ±1  ppbv. 

In  the  case  of  large-scale  model  runs  the  effect  of  moun¬ 
tain  waves  on  the  solid  particle  production  within  the  entire 
northern  polar  vortex  over  the  course  of  the  winter  season 
1999-2000  was  shown  to  be  quite  significant;  Practically  no 
ice  particles  were  seen  as  a  consequence  of  the  synoptic  con¬ 
ditions  alone  whereas  ice  was  observed  for  extended  periods 
of  time  over  the  course  of  the  winter  when  mountain  wave  ef¬ 
fects  were  taken  into  consideration,  regardless  of  the  choice 
of  microphysics.  However,  the  case  where  NAT  nucleation 
above  Tice  was  included  consistently  showed  a  higher  per¬ 
centage  of  trajectories  containing  ice  particles.  When  NAT 
nucleation  above  7)ce  was  included  the  highest  percentage  of 
trajectories  containing  ice  particles  were  around  4%,  whereas 
the  percentage  of  ice  particles  is  the  case  where  NAT  nu¬ 
cleation  above  7)cr  was  not  included  did  not  exceed  0.5%, 
see  Fig.  1 1 .  NAT  particles  were  seen  to  form  readily  as  a 
consequence  of  the  synoptic  conditions  alone,  regardless  of 
the  choice  of  microphysics.  However,  with  the  inclusion  of 
mountain  wave  effects  the  time  period  during  which  NAT 
particles  were  present  became  longer.  As  was  the  case  with 
the  ice  particles,  the  percentage  of  trajectories  containing 
NAT  particles  is  greater  when  NAT  nucleation  above  Tjce 
was  included  (up  to  35%  of  the  trajectories)  than  when  this 
was  not  the  case  (up  to  20%  of  the  trajectories).  Consider¬ 
ing  the  hemispheric  simulations  mountain  wave  effects  seem 
to  have  a  quite  significant  influence  on  the  NAT  formation, 
raising  the  percentage  of  trajectories  containing  NAT  parti¬ 
cles  by  approximately  10%  from  the  warmest  to  the  coldest 
temperature  scenario,  see  Fig.  11.  This  is  contrary  to  what 
was  seen  in  the  two  local-scale  simulations  matching  DC- 


8  flights  examined  here.  A  possible  explanation  could  be 
that  the  synoptic  temperatures  alone  were  low  enough  to  al¬ 
low  for  NAT  formation  during  these  two  flights.  Additional 
local-scale  model  runs  from  relatively  warmer  flight  dates 
(e.g.  early  in  the  winter)  could  help  shed  some  more  light  on 
this. 

Given  the  large  effect  mountain  waves  apparently  had  on 
ice  and  NAT  formation  on  hemispheric  scales  during  the  win¬ 
ter  1999-2000,  the  inclusion  of  mountain  waves  in  models 
seems  important  in  order  to  be  able  to  provide  accurate  esti¬ 
mates  of  the  PSC  load  within  the  vortex.  It  would  be  most 
interesting  to  determine  whether  or  not  the  impact  of  moun¬ 
tain  waves  is  as  pronounced  in  other  winters  as  well  and  such 
an  examination  is  planned  for  a  future  study  along  with  a 
series  of  comparisons  between  hemispheric  model  data  and 
satellite  observations. 

6  Conclusions 

Detailed  lidar  measurements  of  PSCs  made  by  the  NASA 
DC-8  have  been  compared  to  corresponding  model  runs  from 
a  microphysical  PSC  model  using  a  number  of  different  tem¬ 
perature  and  microphysical  scenarios.  It  is  found  that  the  best 
agreement  between  measurements  and  model  data  is  seen  in 
the  case  where  mountain  wave  effects  are  taken  into  consid¬ 
eration  and  where  NAT  nucleation  above  T,ce  is  included  in 
the  microphysics.  Motivated  by  these  findings  the  effect  of 
the  choice  of  microphysics  and  mountain  waves  on  hemi¬ 
spheric  PSC  production  is  examined  in  a  series  of  vortex¬ 
covering  model  runs  extending  over  the  1999-2000  Arctic 
winter  season.  It  is  seen  that  the  localised  mountain  wave- 
induced  temperature  perturbations  have  an  effect  on  the  pro¬ 
duction  of  solid  PSCs  which  extend  beyond  their  local  scale; 
in  the  hemispheric  model  runs  the  percentage  of  trajecto¬ 
ries  containing  solid  PSC  particles  increases  when  mountain 
wave  effects  are  taken  into  account.  Depending  on  the  choice 
of  microphysics,  the  inclusion  of  mountain  waves  is  crucial. 
When  NAT  nucleation  above  Tice  is  not  included  in  the  mi¬ 
crophysics,  practically  no  trajectories  containing  ice  or  NAT 
particles  are  observed  as  a  consequence  of  the  synoptic  con¬ 
ditions  alone,  whereas  both  particle  types  are  seen  to  form 
when  mountain  wave  effects  are  taken  into  consideration. 
When  NAT  nucleation  above  Tu:e  is  included  in  the  micro¬ 
physics,  NAT  particles  can  form  as  a  result  of  the  synoptic 
conditions  alone.  However,  the  amount  of  NAT  and  ice  par¬ 
ticles  increases  as  mountain  wave  effects  are  taken  into  con¬ 
sideration  and  the  percentage  of  trajectories  containing  either 
NAT  or  ice  is  consistently  larger  in  the  scenario  where  NAT 
nucleation  above  T\ce  is  included  in  the  microphysics  com¬ 
pared  to  the  scenario  where  this  is  not  the  case. 
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